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The early Universe may have oscillations, kinations (nonoscillatory evolution), topological defects, 
relativistic and non-relativistic particles at the same time. These components, which may have 
distinguishable density scalings, may decay modulated. In this paper we study generation of the 
curvature perturbations caused by the modulated decay in the multi-component Universe. For our 
analytic calculation we consider the instant-decay approximation. The non-Gaussianity parameter 
is separated into the model-dependent and independent contributions, where the sign of the model- 
independent contribution is determined by the density scalings. We also find that the component, 
whose scaling is very close to the radiation, can generate significant non-Gaussianity from the 
model-independent part, and its sign could be either positive or negative. Importantly, this result 
shows that significant non-Gaussianity can be generated even if the decaying component is already 
dominating the universe. We also consider partial decay, in which a fraction of the component 
decays modulated. 

PACS numbers: 98.80Cq 



I. INTRODUCTION 

Our focus in this paper is the primordial curvature 
perturbation that exists on cosmological scales just 
before they start to enter the horizon. The cosmological 
observations are strongly constraining those cosmological 
perturbations and providing a unique window on the very 
early Universe [l| . 

Generation of the cosmological perturbations begins 
presumably during inflation, when the vacuum fluctua- 
tions of bosonic fields are converted to a classical per- 
turbation. Within this general framework, one can find 
many proposals [ lH 1 0. Il2i. Il3|. 

First recall the SN formalism used to calculate £• To 
define the curvature perturbation £, the energy density 
p is smoothed on a super-horizon scale shorter than any 
scale of interest. Then the local energy continuity equa- 
tion is given by 



dp{x,t) 
dt 



da(x, t) 



a(x,t) dt 



{p(x,t)+p(x,t)), (1) 



where t is time along a comoving thread of spacetime 
and a(t) is the local scale factor. During nearly expo- 
nential inflation, the vacuum fluctuation of each light 
scalar field 4n is converted at horizon exit to a nearly 
Gaussian classical perturbation with spectrum (H/2ir) 2 , 
where H = a(t)/a(t) in the unperturbed Universe. Writ- 
ing the curvature perturbation 



C = S[]n(a(x,t)/a(h)] = SN, 



(2) 



N((f>i(x,U),<f>2(x-,U), 



, t, t*) so that 



C(x,t) = Ni6<f>i(x,t*) + -Nij8^j{x,U) + 



(3) 



where a subscript i denotes d/d<fn evaluated on the un- 
perturbed trajectory. The 8 N- formalism is also valid for 
the evolution after inflation. 

We consider modulation of the decay rate T (<£>), which 
causes inhomogeneous decay of the component (p a ) in 
the multi-component Universe. We are assuming instant 
decay for the calculation @. See FigQ]fbr the basic set- 
ups of the modulated reheating scenario. 

The source of the modulation is the moduli perturba- 
tion caused by an additional light field ip, whose potential 
is assumed to be negligible at the decay. The "seed" per- 
turbation Sip is generated during the primordial inflation. 
At the horizon exit, we consider Gaussian perturbation 
5<p« = ip* — (p*. At the decay, we introduce the function 
<Pd = gi^P*) and the expansion about the Gaussian per- 
turbation Sip = g'Sip* (or equivalently Sip* — Stp/g'). The 
function g explains evolution after the horizon exit [l6| . 

In our scenario of the multi-component Universe, we 
assume that the component p a decays after reheating. 1 
For the first example, we consider simple two-component 



and taking i* to be an epoch during inflation af- 
ter relevant scales leave the horizon, we assume 



1 In the multi-component Universe, the calculation may generically 
depend on isocurvature perturbations between components. Al- 
though we are calculating the perturbations when the isocurva- 
ture perturbations are negligible, our formalism is carefully pre- 
pared so that the mixed perturbations can be calculated within 
the formalism. 
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FIG. 1: Modulation at the transition causes density perturba- 
tions when the decaying component changes its density scal- 
ing. The straight line shows the instant-decay approximation. 



Universe, which has p a oc a 3 (matter) and p r oc a 4 
(radiation) before the decay. 

Later in this paper we are going to extend our ana- 
lytic calculation to the components that may "not" scale 
like matter [13l - ll5| ; typical examples are the cosmological 
defects or the oscillatory / nonoscillatory evolutions Q 
that may be caused by the non-quadratic potentials. 



II. MODULATED DECAY IN THE SIMPLE 
MULTI-COMPONENT UNIVERSE 

First consider the simple (matter + radiation) multi- 
component Universe. The set-ups are similar to the cur- 
vaton. In the curvaton mechanism the significant 
contribution comes from the evolution before the decay; 
while the modulated decay [9( describes the generation 
of the curvature perturbations at the end [l(| [Hj] ■ 

For our analytic calculation, we consider instant-decay 
approximation [9j. The actual transition could be more 
complicated depending on the details of the model pa- 
rameters. For the modulated reheating scenario, the 
idea of the continual decay has been considered by many 
authors [I?], EH- In fact, the Universe in which mat- 
ter is decaying continually into radiation can realize the 
multi-component Universe during the decay. In that way, 
just by shifting the initial conditions, previous numerical 
studies could have been including "the modulated de- 
cay in the multi-component Universe" . In that sense nu- 
merical calculation is versatile. On the other hand, the 
analytic estimation, which is calculated in specific lim- 
its, will help understand the physics behind it. In our 
analytic calculation we found that mo del- independent 
non-Gaussianity is significant both in the single and in 
the multi-component Universe. The mo del- dependent 
non-Gaussianity is calculated from the perturbations of 
r(y?*). Curvature perturbation in models where a light 
scalar field is the curvaton and the moduli that induces 
modulated reheating the same time has been considered 
in Ref.19]. We hope the reader can easily compare our 
calculation with those previous studies. 

In our model, the uniform density hypersurfaces de- 
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FIG. 2: In the left picture we show the densities and their 
scalings in the modulated (r 7^ To) patch. The right picture 
shows the unmodulated (reference) Universe. These pictures 
are sharing Fo (p CT .o and p r .o) and T as the common quantities. 
Due to the different p„-scalings after To, one will find N m 7^ 
and p;,r,m / Pi, re- 



fined for the instant decay is given by 

Pa,T + pr,r = 3M^r 2 



(4) 



where the instant decay occurs at H = T. More specifi- 
cally, in the modulated Universe (see the left picture in 
FigJ2]) we have 



2 r 2 



Pa,r,rn+ Pr,r,m = 3M p T 



(5) 



and in the unmodulated (reference) Universe (see the 
right picture in Figf2]) we have 



2p2 

P 



(6) 



Introducing the mean value To, we define the uniform 
density hypersurface 



Pafi + Prfi = 3M„T, 



2 r 2 
p 1 0> 



(7) 



which is common in the modulated and in the unmod- 
ulated Universe. Without loss of generality, one may 
choose r < Tq for the calculation. 

Using the density scalings, we find in the modulated 
Universe 



Pa,r,m — Pa,0 



a r 



Pr,r,; 



Pr.O 



ar,-. 



a r 



(8) 



which lead to 



Per,0 + Prfi 

Defining local "N m " by 

r*rv 



r 2 

= pr- 0) 



N„ 



H(t)dt, 



(10) 



3 



where H (t) = a/a, one can rewrite the above equation 
as 



Ue-^- + (1 - f a )e-™~ = , (11) 

1 o 

where the coefficient is defined by 

Pa,0 



fa 



Pa.O + Pr,0 



(12) 



In order to compare N m with the unmodulated Uni- 
verse, we find the local equation in the unmodulated Uni- 
verse, 



f a e- iN * + (1 - U)e~ iN * = L. (13) 



Let us include isocurvature perturbations at H — Tq. 
We consider the non- linear formalism [2(|[2l| applied just 
before H = To; 



Cx - SN ini + 



1 



Pa.O 



P<T,0 P° 



4 Jpr.O Pr 



(14) 
(15) 



where p CTj o and p r $ are defined on the uniform density 
hypersurfaces, and pi denotes their mean value. Here 
5N ini denotes the curvature perturbation at that mo- 
ment, which could be significant if the curvaton evolution 
before the decay is significant. 

Finally, the non-linear formalism gives 



Pa,o = Pa,oe 

Pr.O = PaflR 



3(Co— <Wi»i) 



(16) 
(17) 



which leads to 

fa = 



3(CaSN ini ) 



p ffi0e 3(C.-*iV tai ) + p rfie i(CrSN ini ) 

,3(G-*JV ln i) 



Pa,oe 



(18) 



Using the component perturbations C, a and £ r , we find 

r 2 



fad 



3((<,-6N ini -N m ) 



+ (1 - U)e 



4{Cr SN ini -N m ) 



r 2 
1 



1 



where the coefficient is defined by 

Pa,0 



fa 



3M 2 r 2 ' 



(19) 



(20) 



A. First order 

From Eq. (fT9|) . wc find at first order 
ST {1) 



To 



+4(l-/ (T )(cW-«57V i ( V-^i 1) ) (21) 
-9/ CT iVW(C«-^)) 
-16(l-/ ff )iVW(CW-<yJvA)) (22) 



<5r 



(i) 



= U(3&-38NW-4NW) 
+±(l-f a )((W_ SN W- N W) 



-12 faN^(C^-SN^) 

-16(1- U)N$\CP-6NU), (23) 

where the evolution of the perturbations between H = 
Tq and H = T vanishes at this order, because we have 
AKo) =i \r(°) =0 . 



Solving the above equations, we find 



ST 



(i) 



r„ 



N 



(i) 



+(l-r CT )(C«-5<V) 
1<^ 3 



+ (7« ~SN {1} ) 

tA — 2 T 4^ CT ^cr OJV ini ) 



(24) 



(25) 



Pa 



where the coefficients are defined by 

2(At,0 + Pr,o) 

3p CT ,o + 4p r ,o 
3/5 CT ,o + 4/9 r ,o 

Wc thus find the relation 



1 _ 1 

Pa ~ ^ - g^ 



(26) 



(27) 



We are calculating to see the deviation from the un- 
modulated Universe. ^ is not the generation 
of the curvature perturbations in the unmodulated Uni- 
verse. 

Therefore, the curvature perturbation created by the 
modulation (SN^ = N$ - NfP) is calculated as 



6NW = -i 



srtt i <5r 



r 2 r 




(28) 



where the other terms cancel by definition. 

To understand the above result, let us consider the 
perturbation of T. Expanding ifd = gif*), which defines 
if at the decay, we find [22[ 



fd = g + Sif. 



(29) 



4 



where Sip = g'Sip*. In that way the first order perturba- 
tion of the decay rate is calculated as 



dg 



Sip. 



(30) 



g=g 



In the practical calculation ST m and ST^ are identical. 
We thus find 



,v... = ,^ + i)<£, 



(31) 



where ST^ = ST^n = ST^V. In the single-component 
limit (p a — 2/3), we find 



6N™ 



i srw 



6 r 
i r 



6 IV 



Sip, 



(32) 



which reproduces the calculation in Ref. [;9] . 

It is obvious that g is trivial in the slow-roll limit; how- 
ever for more practical estimation one has to calculate 
the model-dependent function g. We are avoiding the 
calculation of <?, since it is highly mode-dependent. 



B. Second order 



Generically, one can expand 



LXJ ^ 



fc=i 



(33) 



where S^ip is a Gaussian random field. In the same way, 
the primordial perturbation can be expanded as 



oo 



fc=2 



(34) 



where £W is Gaussian. Non-linearity parameters are de- 
fined for the adiabatic perturbation £; 

C = C (1) + ^l(C (1) ) 2 + |^l(C (1) ) 3 + -- (35) 

Using the Gaussian quantum fluctuations at the hori- 
zon exit (6<p* ) , we can write [22j | 

ip* — (p* + Sip*, (36) 

which is exact by definition. Again, we write 

Vd = g(ip*) (37) 

and expand it as [22[ 



1. Decay rates 

Before discussing non-Gaussianity of the second order 
perturbations, let us examine the specific expansion of 
the model-dependent decay rates. 



Our first example is 



(39) 



Then, one can expand 

- V°° J-oW ( 
l^k=i k\y \g' <p J 



1 

1+ 2- 



(40) 



We thus find for the expansion T = T a + ST^ 1 ' 



ir( 2 ) 



with the approximation r ~ Tq; 



ST^_ 

r 

ST (2) 



9 x 



1 

Ml 
Ml 



9' 



9" 9 

(g') 2 

g"g 
(g') 2 



(Sip) 2 



(41) 



(42) 



An interesting case would be g < M*, where the 
initial condition is comparable but less than the 
cut-off scale. In that case one can find significant 
fpfL in the conceivable range. Moreover, it is pos- 
sible to find negative contribution from 



rft,) = r(i-i4 

y ^ J 1 2 Ml 



The flip of the sign is very important. 



(43) 



• Second, we consider T oc ip n . The specific form 
beocmes 



^' n M? 
Then, T can be expanded as 



sr^ 
r 

Tn 



Sip 



i i\ , 99 



(<?' 



(44) 

(45) 
(46) 



fini = 9 + Y.T\ 



fc=l 



g 



(fe) f g_^p 
g' v> 



where we wrote g( k ' = d k gjdip\. 



(38) 



Let us summarize the results. Defining 



= A 1 



V r 



(47) 



5 



we find 



.4 



±- 



9 2 



for T(ip) 



g"g 
(g') 2 



1 ± 



1^1 

2 M? 



(48) 



and 



.4 



(n - 1) + 



99 



(g') 2 



for r(v?) 



At 



n Ml 



(49) 



The above calculations are important for the estima- 
tion of the mo del- dependent part of the non-Gaussianity 
parameter Jnl- 



2. f N L 

In order to extract the contributions from the modu- 
lation, we are going to assume ~ C, r ~ 0. We also 
assume 8N- m i ~ for simplicity. 

Then, one can easily expand Eq. (IT91 to find the second 
order perturbations. The expansions used here are 



iN 



and 



l + a(N^ + l -N^ + ...) 
+y(^V (1) + V 2) + ...) 2 + 



'r + ^r( 1 ) + i ( 5r( 2 ) + 



(50) 



(51) 



We find for the second order perturbations 



N. 



(2) 



(O 



r 2 
1 o 



(52) 



Using the relations between the first order perturbations; 



iV« = - 



iV 



(1) 



t 

2p g 
' 1 - 2p 



2 r 
i 

l-2p CT ' 



(53) 



(54) 



and the definition (|47j) . we find 

[^(32-14/,) 4p CT (l + A) 



AT(2) 

ra 
(2) 



iV. 



/// 



(1 - 2p CT ) 2 
4 



(1 - 2p. 
2(1 + A) 



(1 - 2 Va f (1 - 2p CT ) ; 



(55) 



We thus find 

}nl — 



pi (16 - 7/ CT ) - 2p CT - 1 



(1 - 2p CT ) 



5A 



1 



3 l-2p 



(56) 



where the last term depends on A (the model-dependent 
contribution). In the single-component limit (/ CT — > 1 
and p a — > 2/3), we find a simple formula 



fNL = 5-5A 



(57) 



Note that the model-independent contribution /jvx = 5 
is showing an interesting result. Note also that r oc tp 3 
gives A ~ 1 when the evolution of </? is trivial, which may 
lead to the cancellation (fNL — 0). 

In the opposite limit, / CT — >• leads to p CT — » 1/2. In 
that limit we find 



fNL oc 



1 



(1 - 2p CT )- 



> 1. 



(58) 



III. 



HIGHER POTENTIAL OR TOPOLOGICAL 
DEFECTS 



For the scalar potential of the form V(a) oc er", the 
energy density of the scalar-field oscillations decreases as 
P(j oc ci n + 2 when the oscillations are rapid compared with 
the expansion rate [23]. Alternatively, one may choose 
topological defects for the decaying component, which 
may scale like p a oc a*. NO (Non Oscillatory) motion can 
lead to a different density scaling [1J|. Here the scaling 
is approximately defined at the time of the decay; there 
is no need to find exact scale-dependence that is valid 
during the whole evolution. This point might be crucial 
for the practical investigation. 

For our purpose, we consider the component that scales 
like p a oc a~( 4+e "). Here e n = —1 corresponds to the 
sinusoidal oscillation for the quadratic potential, whose 
energy density scales like p a oc a~ 3 . Note that e„ > 
is not excluded in our calculation; we will show that this 
may change the sign of /jvl- 

In order to include the isocurvature perturbation at 
H = Tq, we consider the component perturbations de- 
fined by 



(a = SN ini + 



Cr = 5N ini + 



1 



4 + e r , 



Pa.O P° 

! r Pr -° Clp r 



^ (59) 



Pr,0 



Pr 



(60) 
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Then we find 

J CTe (4+^)(C„-'5JV 1 » 1 -A r m ) + (1 _ f^ e MCrSN ini -N m ) = £^ 

^0 

^ e (4+e„)(f <r -*JV iDi )-4JV^ + / 1 _^ e 4(C r -«JV tal -iV^) = £_ 

^0 

(61) 

Again, the isocurvature perturbations and <5Aj n i does not 
contribute to the modulation at first order. They will be 
neglected in the second order calculations. 



A. First order 
From Eq. (|61[) . we find at first order 

<5r (1) 



To 



(4 + e n )/ CT A«-4(l-/ ff )AW (62) 



- -^N^-4(l-f a )N£\ (63) 



To 

Solving the above equations, we find 

<5r (1) 



To 

^(l) = 1 61 

* 2 r 
where the coefficient is defined by 

2(/Oa,o + /5r,o) 



Per,, 



(4 + e n )p ff) o + 4p rj0 



(64) 
(65) 

(66) 



Therefore, the curvature perturbation created by the 
modulation is given by 



6NW ee A« atW 



2/ r n 



(67) 



Obviously, generation of the curvature perturbation is 
possible when e n ^ (i.e, when two components (p a and 
p r ) are distinguishable in their scaling relations). 



B. f N L 

Again, we find for the second order perturbations 
'(4 + e„) 2 /. - ^ ' <^ 2 



Pa 



Pa, 7i 



A 



(2) 



r/i 



•8(1-/,)) (a«) 



lA a (<5r( 1 )) 2 + r ff( 2 



A. 



(68) 



Using the relations 

At(D 

r(l) 



2p CT ,r 



^ = - 



1 - 2p CTj „ 

i <rrw 

2 r 

1 



-<5A« 



l-2 Po 

and the definition ([4"T[). we find 



-<5A«, 



(69) 



(70) 



Ar(2) = 



{2(4 + c„) 2 / ff + 32(1- /„)}- 



Pa,n 



(1 - 2p ff ,„) 5 



4 P(T , ra (l + A) 
' (1 - 2?w) 2 
4 

(1 - 2p CT ,„) 2 " (l-2 Ptr , T 



2(1 + A) 



<5A« 



(71) 



We thus find 



INL 



r p a (4 + ej 2 /,, + 16^(1 - - 2p CT - 1 



(1 - 2p ff ) 2 



5.4 



1 



3 l-2 Pa ' 



(72) 



where the last term gives the model-dependent contribu- 
tion. The single-component limit is given by f a — > 1 
and p a — > 2/(4 + e n ), where we find significant non- 
Gaussianity 



Jnl 



5(4- 



5A, 



(73) 



which shows that the sign of the first term (model- 
independent contribution) is determined by e n - We find 
positive sign for e n < 0, while it goes negative when 
e n > 0. Moreover, the result shows significant non- 
Gaussianity in the model-independent part, which could 
be very large even if the decaying component is dominat- 
ing the Universe. 

In the opposite limit, f„ — > and p a — > 1/2, we find 



Snl oc 



1 



;i - 2 Pff ) 2 



» i, 



(74) 



as expected. 



IV. PARTIAL DECAY 



More practically, there could be a moment when a frac- 
tion of the matter component decays modulated and the 
decaying component does not have significant interac- 
tion with the other (matter) component. This could be 
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realized when the non-relativistic matter contains parti- 
cles that belong to the hidden sector. 

For the multi-component Universe that contains both 
matter and radiation, the uniform density hypersurfaccs 
defined for the partial decay is given by 



2^2 



Pa,T + &po,T + Pr,T = SM'T 



(75) 



where the decaying component is denoted by Ap CT .r- 
Ignoring isocurvature perturbations and SNim, we find 



fae 



fae 



-3N„ 



+ /Ae 



-3N„ 



+ 0--U- U)e 



-4jV„ 



+ /Ae- 4 ^ + (1 - U - /a)< 



-4iV„ 



where the coefficients are defined by 

Po-,0 



U 

Fa 



Pafi + Ap a0 + Pr,0 
Pc.Q + Ap o + PrM 



First order 



We find at first order 



To 



= -3/ ff iV« - 3/aA^ 
-4(1 - U - h)N$ 

i o 

-4(1 f(j /aX 1 ) 
Solving the above equations, we find 

5TW 



7V« = 



iVi 15 = 



-PA- 



r 



* r 

where the coefficients are defined by 

2(p <T) o + Ap CT . + p r , ) 



PA 



3/Oa,0 + 3A~P<t,0 + 4 Pr-,0 
2(p<7,0 + APtrfl + Prfi) 

3p<j : o + 4Ap + 4p r>0 



r 2 
1 {) 



r 2 ' 

1 



(76) 

(77) 
(78) 



(79) 

(80) 

(81) 
(82) 

(83) 
(84) 



Therefore, the curvature perturbation created by the 
modulation is 



<WV« = atM _ N 



(1) 



(-PA+ £>a) 



5TV 
To 



-PATA- 



(85) 



where the last approximation is valid when rA *C 1. Here 
the coefficient is defined by 



Ap a fi 



3p afi + 3Ap o + 4p rj0 



B. Second order 



(86) 



We find for the second order perturbations 



pa (z-l(U + fS) (n£> 



-PA- 



r 2 

1 



(87) 



£>A f 8 - -fa) 

(<5r«) 2 + rvr( 2 ) 



Using the relations 

±y m 

r(l) 



r 2 

1 



-PA- 



r n 



PA 



PA - i?A 
A-fT~ 



PA - i> A 

and the definition (|4"T|) . we find 

Pi(16 - 7/ CT - 7/ A ) Pa(1 + ^) 



(89) 



(90) 



tv( 2 ) 

m 
(2) 



2(PA " t>A? 



(pa - ^a) 2 



A'. 



j> A (16-7/ CT ) 

2(PA - i>A? (PA - t> A ) 2 



(SN^Y . (91) 



Then, /tvl is calculated from 



5 a4 2) -< } 



(92) 



V. CONCLUSION AND DISCUSSION 

The early Universe after inflation may have many com- 
ponents. They could be oscillations, topological defects, 
relativistic and non-relativistic particles. If those compo- 
nents are decaying into radiation, there could be gener- 
ation of the curvature perturbations. In this paper, the 
mechanism of the modulated decay has been considered 
for the multi-component Universe. 

In this paper we found the basic formulation, which 
is useful in calculating modulated decays in the multi- 
component Universe. 
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We found useful results, in which the non-Gaussianity 
parameter is separated into model-dependent and inde- 
pendent terms. Interestingly, the sign of the model- 
independent Jnl may flip depending on the scaling of 
the decaying component. 

We find that the component whose scaling is very close 
to the radiation (e„ ~ 0) could generate significant non- 
Gaussianity in the single-component limit. Importantly, 
the sign could be either positive or negative. In that way, 
the modulated decay of the oscillations may crucially de- 
pend on the amplitude at the decay. For instance, con- 
sider the potential for the oscillations given by 



\rl \ ~ 1 2 2 , ^4 4 A 6 if 

V(a)^-ma + -a + 



(93) 



If the oscillations decay when ip e is dominant, one will 
find Jnl < 0. If the oscillations decay when (p 4 is dom- 
inant, one will find \$nl\ 1, where the sign could be 
either positive or negative. The scaling of the density 
changes during the oscillations. One will find conven- 
tional result when the quadratic term is dominating. In 
the intermediate region one may find the density scal- 
ing p a oc a _fc , where (effectively) 3 < k < 6 is possible. 
As the result, in the practical calculations the curvature 
perturbations and the non-Gaussianity may depend cru- 
cially on the amplitude of the oscillations. 
Note added: While finalizing this paper, we found a 
couple of papers [24j ] which has some overlaps with our 
models. 
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Appendix A: Non-Linear formalism and the 
curvaton mechanism 



In this section we first review the basics of the non- 
linear formalism and then compare our results with 
Ref . [2~i| . Notations in Ref.[24| are discriminated by the 
subscripts "LT-A". 

The non-linear formalism in the curvaton mechanism 
is given by the formula 



Ca = 5N + 1 -ln(^ 

3 \pcr 



(Al) 



Cr 



SN +]ln * 

4 \Pr 



(A2) 



Here SN is the perturbation of N, which is measured 
between two hypersurfaces. Usually the one is the flat 
hypersurface and the other is a uniform density hyper- 
surface. Besides SN, we have to define the other quan- 
tities (p a , p r ) and (p a , p r ). Following the definition of 
SN, those quantities are defined on the uniform density 
hypersurface on which SN is defined. 

We thus find at the uniform density hypersurface H = 

Ua = SN A (t A ) + \ In (PzA^Al) , (A3) 
3 V p„, A (t A ) J 

Cr,A = 5N A (t A ) + \lJ P - A{ f A A . (A4) 
Solving these equations we find 



Pa, A 



Pa,A<i 



Pr,A — Pr,Ae 

The trivial identity is 

P<r,A + Pr,A 
Po,A + Pr,A 



4(Cr,ASN A ) 



= 1, 



(A5) 



(A6) 



where p a ,A and p r>A can be replaced using Eq. (|A5[) . We 
find the equation 

f a , A e 3 ^- m ^ + (l-f„, A )e^- SN ^ = 1,(A7) 

where the ratio is defined by 

Pa, A 



fa,A = 



Pa,A + Pr,A 



(A8) 



We find at first order 



SN A = r CT:j4 Ccr,A + (1 - »V,o)Cr,o 

= SN A 

3 V Pa. A / 4 



Pr,A 



(A9) 



The trivial identity is 



rA In ( ^ 

Pa, A 



4 \Pr,A, 



For the expansion Spi = p t — pi, the above equation gives 
the obvious identity 



Spa,A + Sp r , A = 0. 



(All) 



One may evaluate the non-linear formalism away from 
H = H A . (See Figl3l) Choosing another hypersurface 
at H — Hb, one can evaluate 

SN B = r a , B (a,B + (1 - r a . B )Cr,B 

= r a . B Ca, A + (1 - r atB )Cr, A , (A12) 
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(P 



SN 



A 



Pr,A); (pa,Ai Pr,A) 



SN 



curv 



H = H 



A 



SN 



B 



H = H 



B 



(pa,B, Pr,B)] (Pa,B, Pt,b) 



FIG. 3: Component perturbations are denned on the uniform 

density hypersurfaces. For instance, C,a,A = SNa + \ In J 

is defined using SNa (measured from the flat hypersurfaces 
to H = Ha), pa, a and its mean value pa, a (both are de- 
fined on H = Ha)- Note that the isocurvature inhomogeneity 
Spa,A(x, tA) = Pa,A(x,tA)~ Pa,A(tA) is possible in the uniform 
density hypersurfaces when the Universe is multi-component. 



where the constancy of the component perturbations 
Cm = d,B has been used. 

Note that 5N CUIV = SNb — SNa gives the "evolution" 
between the two hypersurfaces. We thus find for r a _B 3> 
TV, a: 



5N C , 



{r<r,B - r a ,A)(cr,A 
Spa,A 



r a ,B 



{r a ,B - r a> A)C,r,A 

(A13) 



Note that SNa does not appear in the final result because 
of the obvious cancellation. 

If one uses Ha at the beginning of the curvaton oscil- 
lation and Hb at the decay, SN CUTV gives the evolution of 
the curvature perturbation in the conventional curvaton 
mechanism. 

The curvature perturbation generated before the cur- 
vaton oscillation can be included as 6N{ n f ~ SNa- Then 
one will find the post-decay curvature perturbation 



SN E 



SN, 



SN C 



(A14) 



Why difficult? 



The formalism that can be applied for the modula- 
tion at the end has already been discussed by Enomoto- 
Kohri-Matsuda(EKM) in Ref.[ll|. When the decay is 
modulated, the non-linear formalism at the end can be 
given by 



SN E 



111 



Po_ 

Po 



ilnf^ 

4 \po,d 



= SN E 

Cr,E = SN E 

where we defined 



In 



P<T,E 

Po 
Pr.E 
Pr 



Po,d 
Po,E 



(A15) 



(A16) 



Here the subscript "E" means that the quantities are 
evaluated at t = t E , which denotes the time after the 
decay. The quantities denoted by "bar" is placed before 
the decay. 

Again, the trivial identity 



Po.E + Pr.E 
Per + Pr 



Hi 

H 2 



gives 



3(G 



-6AT E -A E ) 



+ (1 - U)e 



(A17) 



Hi 
(A18) 



It is possible to identify He = T and H = Tq to find 
Ae = 0, however in that case the perturbation of He = T 
is a problem. We thus need some trick to find the exact 
result. 

At this moment we have two solutions for the problem. 
One is discussed in this paper, and the other is discussed 
in Ref. 24] by two groups. For instance, Langlois and 
Takahashi introduced a new parameter SNb to define 



C, = SN D + hn(^ 

3 V p a 



(A19) 



where SNb cannot be identified with the curvature per- 
turbation while Co- *s identical to the conventional com- 
ponent pertubation. Below, we will take a closer look at 
these definitions. 



2. How to understand the calculation of Ref. [2 
terms of our set-ups 

We are going to show the obvious correspondences be- 
tween the quantities defined in Ref. 2M ano - our set-ups 
presented in FigJ^J 

In Ref. [24| , they have defined the non- linear formalism 



( a = SN D + i In 



Pojtp) 
Pa{t D ) 



(A20) 



where the uniform density hypersurfaces defined by to 
(H = r) and Ib {H = f ) are not identical. SNb is 
defined for the hypersurface H — T, but it is not the 
usual curvature perturbation. (See Fig]?]) 
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Unmodulated 
Decay 




h = r 



-iln(l + <5 r 



^L=r 

Pa (to) Modulated 
Decay 

FIG. 4: Unlike Fig[3] component perturbations of ref.[24| 
are denned using multiple hypersurfaces. For instance, 

CxI(lt-a) = SN D + | In ( ^(fp) ) is denned using both 
H = F and if = T. As a result, -|ln(l + 5r) has to be 
subtracted from 8Nd to get the curvature perturbation £. 



In our formalism, these quantities are corresponding to 

8N D \ {LT _ A) O SNi n i + N m 

Pa(tD)\(LT-A) ±* PaT,m — Pafi£-~ 3Nm 
P<r(tD)\{LT-A) *+ P<r,0- (A21) 

We thus find 

Crl(LT-A) ^ Ct- (A22) 

They also defined 

(1 + S r ) 2 =T 2 /f 2 , (A23) 

which gives the correspondence 

(1 + <5 r ) 2 1 { lt~a) oe~ 4J V. (A24) 

Finally, they have defined the post-decay curvature per- 
turbation 



( = 5N D + ~hi(l + 5 r ), 



(A25) 



which corresponds to 

C\(lt-a) SN ini + N m - N+. (A26) 

In our calculation the curvature perturbation generated 
by the modulated decay is given by 



Therefore, the correspondence is obvious between our cal- 
culation and Ref . [24l] . 

However, in finding the curvaton contribution they 
evaluated 



r r 



(A28) 



where the first and the last terms are originally given by 



Cr + = r( a + (1 - r)Cr = SN ini , 



(A29) 



where (5iV; n j is defined in this paper. If the curvaton hy- 
pothesis is valid and the component perturbations are 
constant, one may evaluate the component perturbations 
at H = H osc : 

(a = Utosc) = SN ini + \ In (A30) 

*J Pa t,^osc ) 

where 8N- m f denotes the curvature perturbation just at 
the beginning of the oscillation. Substituting the com- 
ponent perturbations (defined at H osc ) into the above 
equation, one finds 



Cr + g<5 = r( a (t osc ) + (1 - r)( r (t osc ) 

c Ar i r i Pa(t osc ) 
= 6Ni ni + - In ■ 



3 Pa(iosc) 

l-r p r (t osc ) 
— a — ^7 V ' 

4 Pr \Lqsc) 



(A31) 



Although trivial deformation is needed, it is easy to find 
that the result is consistent with Ea. llAlty without using 
approximations. 

Using the above formula, they started perturbation 
with regard to the perturbation of S. For instance, Lan- 
glois and Takahashi considered for the "curvaton pertur- 
bation" 

S = *(Cr-Cr) = ^-K + lK> (A32) 

a o~ 6 a 6 
and for the "inflaton perturbation" 



C r = -rty ^ 6N in{ . 



(A33) 



<W mod = N m - N+. 



(A27) 



These definitions (in fact they are approximations valid 
only when one can disregard Sp r /p r ) are based on the 
usual curvaton hypothesis. 

Although in this paper we were not focusing on the 
mixed perturbations, our formalism is carefully prepared 
to explain the modulated decay in the multi-component 
Universe, where the curvaton hypothesis could not apply. 



11 



[1] D. H. Lyth, A. R. Liddle, "The primordial density pertur- 
bation: Cosmology, inflation and the origin of structure," 
Cambridge, UK: Cambridge Univ. Pr. (2009) 497 p. 

[2] D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, 
"A New approach to the evolution of cosmological per- 
turbations on large scales," Phys. Rev. D 62, 043527 
(2000) [astro-ph/0003278| ; T. Matsuda, "Delta-N formal- 
ism for the evolution of the curvature perturbations in 
generalized multi-field inflation," Phys. Lett. B 682, 163 
(2009) [arXiv:0906.2525l [hep-th]]: T. Matsuda, "Elliptic 
Inflation: Generating the curvature perturbation without 
slow-roll," JCAP 0609, 003 (2006) [hep-ph/0606"i"37] . 

[3] F. Bernardeau and J. -P. Uzan, "NonGaussianity in 
multifield inflatio n," Phys. Rev. D 66, 103506 (2002) 
[hep-ph/0207295| ; L. E. Allen, S. Gupta and D. Wands, 
"Non-gaussian perturbations from multi-field inflation," 
JCAP 0601, 006 (2006) [astro-ph/0509719| ; H. As- 
sadullahi, H. Firouzjahi, M. H. Namjoo and D. Wands, 
"Curvaton and the inhomogeneous end of inflation," 
larXiv:1207.7006l [astro-ph.CO]. 

[4] T. Matsuda, "Modulat ed Inflation," Phys. Lett. B 665, 
338 (2008) |arXiv:0801.2648l [hep-ph]]; T. Matsuda, 
"Modulated inflation from kinetic term," JCAP 0805, 
022 (2008) [arXiv:0804.3268l [hep-th]] . 

[5] T. Moroi, T. Takahashi, "Effects of cosmological mod- 
uli fields on cosmic microwave background," Phys. Lett. 
B522, 215-221 (2001). [hep-ph/01 10096] , D. H. Lyth, 
D. Wands, "Generating the curvature perturbation 
without an inflat on," Phys. Lett. B524, 5-14 (2002). 
[hep-ph/01 10002] . 

[6] K. Dimopoulos, K. Kohri and T. Matsuda, "The 
hybrid curvaton, " Phys. Rev. D 85, 123541 (2012) 
|arXiv: 1201.60371 [hep-ph]] : T. Matsuda, "Free light fields 
can change the predictions of hybrid inflation," JCAP 
1204, 020 (2012) |arXiv:1204.0303l [hep-ph]] . 

[7] K. Dimopoulos, K. Kohri, D. H. Lyth and T. Mat- 
suda, "The inflat ing curvaton," JCAP 1203, 022 (2012) 
|arXiv:1110.295T1 [astro-ph.CO]] ; 

[8] K. Kohri, C. -M. L in and T. Matsu da, "PBH from the 
inflating curvaton," arXiv:1211.2371 [hep-ph]. 

[9] G. Dvali, A. Gruzinov and M. Zaldarriaga, "A 
new mechanism for generating density perturbations 
from inflation," Phys Rev. D 69, 023505 (2004) 
[arXiv: a stro-ph/030359 1] . 
[10] D. H. Lyth, "Generating the curvature perturbation 
at the end of i nflation," JCAP 0511, 006 (2005) 
astro-ph/0510443 ; F. Bernardeau, L. Kofman and J. - 
P. Uzan, "Modulated fluctuations from hybrid inflation," 
Phys. Rev. D 70, 083004 (2004) [astro-ph/0403315| . 
[11] S. Enomoto, K. Kohri and T. Matsuda, "Non- 
Gaussianity in the unified curvaton mechanism : The 
generalized curvaton mechanism that comprehends mod- 
ulation at the transition," arXiv:1210.7118 [hep-ph]. 



[12] T. Matsuda, "Cosmological perturbations from inho- 
mogeneous preheating and multi-field trapping," JHEP 

0707, 035 (2007) [arXiv:0707.0543l [hep-th]]. 

[13] T. Matsuda, "Cosmological perturbations from an inho- 
mogeneous phase transition," Class. Quant. Grav. 26, 
145011 (2009). |arXiv:0902.4283l [hep-ph]]: M. Kawasaki, 
T. Takahashi and S. Yokoyama, "Density Fluctuations 
in Thermal Inflation and Non-Gaussianity," JCAP 0912, 
012 (2009) [arXiv:0910.3053l [hep-th]]. 

[14] G. N. Felder, L. K ofman and A. D Linde, Phys. Rev. D 
60, 103505 (1999) [hep-ph/9903350| ; T. Matsuda, JCAP 

0708, 003 (2007) [arXiv:0707.1948T [hep-ph]]. 

[15] T. Matsuda, "Topological curvatons," Phys. Rev. D 72, 
123508 (2005) [hep-ph/0509"063] . 

[16] D. H. Lyth, "Can the curvaton paradigm accommodate 
a low inflation scale?," Phys. Lett. B 579 (2004) 239 
[hep-th/0308110] . 

[17] K. A. Malik and D. H. Lyth, "A numerical study of non- 
gaussianity in the curvaton scenario," JCAP 0609, 008 
(2006) [astro-ph/0604387| . 

[18] M. Zaldarriaga, "Non-Gaussianities in models with 
a varying inflaton decay rate," Phys. Rev. D 69, 
043508 (2004) [astro-ph/03"0 6006 ; T. Suyama and 
M. Yamaguchi, "Non-Gaussianity in the modulated re- 
heating sce nario," Phys. Rev. D 77 (2008) 023505 
[arXiv:0709.2545l [astro-ph]]; K. Kohri, D. H. Lyth 
and C. A. Valenzuela- Toledo, "Preheating and the 
non-gaussianity of the curvature perturbation," JCAP 
1002, 023 (2010 ) [Erratum-ibid. 1009, E01 (2011)] 
[arXiv:0904.0793l [hep-ph]]; M. Cicoli, G. Tasinato, 
I. Zavala, C. P. Burgess and F. Quevedo, "Modulated 
Reheating and Large Non-Gaussianity in String Cosmol- 
ogy," JCAP 1205, 039 (2012) |arXiv:1202.4580l [hep-th]]. 

[19] K. -Y. Choi and O. Seto, "Modulated reheating by cur- 
vaton," I arXiv: 1204.14191 [astro-ph.CO]. 

[20] D. H. Lyth, K. A. Malik and M. Sasaki, "A General proof 
of the conservation of the curvature perturbation," JCAP 
0505, 004 (2005) [astro-ph/04lT220] ; 

[21] D. Langlois, F. Vernizzi and D. Wands, "Non-linear 
isocurvature perturbations and non-Gaussianities," 
JCAP 0812, 004 (2008) |arXiv:0809.4646l [astro-ph]] . 

[22] M. Sasaki, J. Valiviita and D. Wands, "Non-Gaussianity 
of the primordial perturbation in the curvaton model," 
Phys. Rev. D 74, 103003 (2006) [astro-ph/0607627| . 

[23] M. S. Turner, "Coherent Scalar Field Oscillations in an 
Expanding Universe," Phys. Rev. D 28, 1243 (1983). 

[24] D. Langlois and T. Takahashi, "Density Pertur- 
bations from Modulated Decay of the Curvaton," 
larXiv: 1301.33191 [astro-ph.CO]; H. Assadullahi, 
H. Firouzjahi, M. H. Namjoo and D. Wands, "Modulated 
curvaton decay," arXiv:1301.3439 [hep-th]. 



